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ABSTRACT 

We propose that the gravitational collapse of supermassive objects 
(M 10^ M0), either as relativistic star clusters or as single supermassive 
stars (which may result from stellar mergers in dense star clusters), could 
be a cosmological source of 7-ray bursts. These events could provide the 
seeds of the supermassive black holes observed at the center of many galaxies. 
Collapsing supermassive objects will release a fraction of their huge gravitational 
binding energy as thermal neutrino pairs. We show that the accompanying 
neutrino/antineutrino annihilation-induced heating could drive electron/positron 
"fireball" formation, relativistic expansion, and associated 7-ray emission. The 
major advantage of this model is its energetics: supermassive object collapses are 
far more energetic than solar mass-scale compact object mergers; therefore, the 
conversion of gravitational energy to fireball kinetic energy in the supermassive 
object scenario need not be highly efficient, nor is it necessary to invoke 
directional beaming. The major weakness of this model is difficulty in avoiding 
a baryon loading problem for one dimensional collapse scenarios. 

Subject headings: gamma rays: bursts - cosmology: observations and theory 
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1. Introduction 

In this letter we propose that the collapse of supermassive objects and the associated 
neutrino/antineutrino annihilation could give rise to high redshift (cosmological) 7-ray 
bursts (GRBs). This model could alleviate vexing problems associated with the energetics 
of conventional stellar remnant-based scenarios. We define a supermassive object to be a 
star or star cluster that suffers the general relativistic Feynman-Chandrasekhar instability 
during its evolution. This corresponds to objects with initial masses M ^ 10'^ Mq, i.e., 
those which may leave black hole remnants with masses M ^10^ Mq. 

Detections of absorption and emission features at a redshift z = 0.835 in the spectral 
observation of the afterglow of 7-ray burst GRB970508 (Metzger et al. 1997a,b) and at 
redshift z = 3.42 in the host galaxy of GRB971214 (S. Kulkarni et al. 1998) have established 
that at least some of the GRB sources he at cosmological distances. Observations show 
that the total energy in gamma rays associated with a GRB at cosmological distances 
is ~ 10^^ erg to ~ 10^^ erg when a An sohd angle coverage is assumed (Fenimore et al. 
1993; Wijers et al. 1997; Kulkarni et al. 1998). Catastrophic collapse events, such as 
neutron-star/neutron-star mergers (Paczyhski 1986; Goodman 1986; Eichler et al. 1989), 
neutron-star/black-hole mergers (Mochkovitch et al. 1993), failed supernovae (Wooseley 
1993), "hypernovae" (Paczyhski 1997), collapse of Chandrasekhar-mass white dwarfs (Usov 
1992), have been touted as natural candidates for cosmological GRB sources. Fireballs 
created in these collapse events could accelerate material to the ultra-relativistic regime, 
with Lorentz factors P = Ee/nieC^ ^ 10^ (Paczyhski 1986, Goodman 1986, Rees & 
Meszaros 1992, Meszaros & Rees 1992). The kinetic energy in these fireballs could then 
be converted to 7-rays possibly via the cyclotron radiation and/or the inverse Compton 
processes associated with ultrarelativistic electrons. In these models, the energy loss of the 
shock(s) propelled by the fireball would produce the afterglow associated with a GRB event 
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(Waxman 1997). 

There are, however, problems for these stellar remnant-based models if the GRBs 
originate from high redshift events. The total gravitational binding energy released when a 
~ 1 M0 configuration collapses to a black hole (or into a pre-existing larger black hole) is 
only ~ 10^^ erg. Calculations have shown that it is very difficult to power a GRB of energy 
~ 10^^ erg (Wijers et al. 1997), or an afterglow with a similar energy (Waxman 1997; Dar 
1997) with such a collapse scenario, unless the 7-ray emission and the blast wave causing 
the afterglow are highly coUimated (improbably highly coUimated in the case of very high 
redshift events). 

This energetics problem can be avoided in the supermassive object collapse model 
suggested here. Collapse of such large mass scale objects could result in prodigious 
gravitational binding energy release. Some of this gravitational energy is radiated as 
thermal neutrino/antineutrino pairs (Fuller, Woosley, & Weaver 1986, hereafter FWW; 
Fuller & Shi 1997) whose annihilations into electron/positron pairs could create a fireball 
above the core that generates 7-rays. There is no direct evidence for supermassive stars 
ever having been extant in the universe. However, it has been argued that their formation 
could be an inevitable result of the collapse of ~ lO^Mo to lO^Mo primordial clouds (the 
baryon Jean's mass at early epochs, see Peebles & Dicke 1968, and Tegmark et al. 1997) at 
high redshifts in which cooling was not as efficient as in clouds contaminated with metals, 
or more likely, as a result of stellar mergers associated with ^ 10^-IO^Mq relativistic star 
cluster collapse (Hoyle & Fowler 1963; Begelman & Rees 1978; Bond, Arnett, & Carr 1984; 
FWW; McLaughlin & Fuller 1996). The flow chart for supermassive black hole production 
suggested by Begelman & Rees (1978) includes several pathways whereby supermassive 
stars are formed in the central region of the collapsing cluster. Further, supermassive black 
holes apparently are ubiquitous in the universe. They are invoked as the central engines of 
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Active Galactic Nuclei (AGNs) and quasars, and are inferred to be in the centers of nearby 
galaxies (van der Marel et al. 1997). 

We note that Prilutski and Usov (1975) have previously tied GRBs to magneto-energy 
transfer during collapses of supermassive rotators (~ IO^Mq) postulated to power AGNs 
and quasars. Here we propose a different energy transfer mechanism (neutrinos) based on 
objects not necessarily tied to AGNs or quasars, but which could possibly be related to the 
birth of the supermassive black holes that power them. 



2. Fireballs from Supermassive Object Collapse 

Supermassive stars will suffer the General Relativistic (Feynman-Chandrasekhar) 
instability, either at or before the onset of hydrogen burning (c.f., FWW) in the case of 
quasi-statically contracting objects, or immediately upon formation as in the case where 
stellar mergers produce them. As such a star collapses, the entropy per baryon is slightly 
increased by nuclear burning, but then is reduced by neutrino pair emission. Though 
initially the whole star can collapse homologously, as the entropy is reduced only an inner 
"homologous core" can continue to collapse homologously (FWW). It is this homologous 
core that will phmge through an event horizon as a unit to make a black hole. The 
mass of the homologous core, M^^ = M^^ / 10^ Mq, can be much smaller (possibly by an 
order of magnitude or more) than the mass of the initial hydrostatic supermassive star. 

The collapse to a black hole of a supermassive star with a homologous core mass 
M^*-' will have a characteristic (prompt) Newtonian gravitational binding energy release of 
~ £"3 ~ 10^^ M^*^ erg. During the collapse, neutrino emission will ensue from e^'^-annihilation 
in the core. The emissivity of this process scales as the core temperature to the ninth power 
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(Dicus 1972). As a result, most of the gravitational binding energy removed by neutrinos 
will be emitted very near the point where the core becomes a black hole, and on a timescale 
characterized by the free fall time (or light crossing time) of the homologous core near 
the black hole formation point. We employ a characteristic free fall collapse time scale of 
ts ~ M|^*-' sec, and a characteristic radius (the Schwarzschild radius) of Tg fa 3 x 10^°M|^'-' cm. 
For a core mass ^ 10^ Mq the neutrinos will not be trapped in the core and will freely 
stream out. For a smaller core mass, the neutrino diffusion timescale will be long compared 
to the free fall timescale and so neutrinos will be trapped in the core. Neutrino emission in 
this latter case will be from a "neutrino sphere" at the edge of the homologous core. 

In general it is difficult to estimate the range of initial stellar masses that will give 
rise to a given range of homologous core masses, though there is a clear hierarchy at each 
evolutionary stage. We therefore guess that the initial star cluster masses will be in the 
range 10^ Mq to 10^ Mq, while the subsequently produced supermassive stars will have 
masses M™* fa 0.1 to ~ 1000, while the corresponding homologous core masses will he in 
the range M|^^ fa 10"^ to ~ 10. Figure 1 shows a flow chart for the collapse of supermassive 
objects. 

The neutrino luminosity can be crudely estimated from the product of the neutrino 
energy emissivity (Schinder et al. 1987; Itoh et al. 1989) near the black hole formation point 
and the volume inside the Schwarzschild radius, i.e., 4 x 10^^ (Tg^'^'^^)^ (47rrg/3) erg/sec. 
Here Tg'^^'^ is the characteristic average core temperature near the black hole formation 
point in units of 10^ K. For a spherical non-rotating supermassive star we can show that 



where aschw is the ratio of the final entropy per baryon to the value of this quantity in the 
initial pre-coUapse hydrostatic configuration, and gg ~ gb + 7/8gf ~ 11/2 is the statistical 
weight of relativistic particles in the core. Since for spherical non-rotating supermassive 
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stars M^^'^/Mf^ ^ V5.5/2agJ^^ (FWW), we can conclude that Tg^^^^ ^ l?>{Mf^)-^/'^ . The 
characteristic neutrino luminosity is then 

L^p ~ 4 X 10^^ (Tf i^w)9 (47rr3/3) erg/sec ^ 5 x IQ^'^ {Mf^)-^/'^ erg/sec. (2) 

Since 70% of the neutrino emission is in the v^.v^ channel, the characteristic luminosity 
of Ve or z/e is Ly^ = Lp^ ^ 0.35L,yp. (This estimate of L^p is a factor of ~ 10 above an 
appropriately scaled version of the Woosley, Wilson and Mayle (1986) result for a M™* = 5 
configuration; part of the difference is attributable to the employment of different neutrino 
emissivities, and the remainder may result from different core temperatures.) 

The copious uu emission during the collapse can create a fireball above the homologous 
core by — > e~^e~ . Clearly, the neutrino luminosities will suffer gravitational redshift which 
will degrade the total energy deposition above the star, though this will be compensated by 
increased i/z/-annihilation from gravitational bending of null trajectories (Cardall & Fuller 
1997). A detailed calculation of these two effects is beyond the scope of this paper, but 
we do not expect the combination of them to change our order-of-magnitude estimates 
significantly. The energy deposition rate per unit volume from the annihilation at a 
radius r above a spherical shell of thermal neutrino emission with a radius R^, is then 
(Goodman, Dar, & Nussinov 1987; Cooperstein, van den Horn, & Baron 1987) 

^""'^'^ - W Rt " ' \-W) WV ^ ^ 

Here Gp is the Fermi constant, L is the luminosity of the neutrinos/antineutrinos, and 

the brackets denote averages of neutrino energy or squared-energy over the appropriate 

neutrino or antineutrino energy spectra (see Shi & Fuller 1998). The phase space and spin 

factors are K 0.124 (0.027) for v = v^, (z/^,//,-), and the radial dependence of the energy 

deposition rate is $(x) = (1 — xf [x^ + 4x + 5), with x = [1 — (-R;,/r)^]^/^. 

The characteristic neutrino luminosity Lyy in eq. (Q) could be an underestimate of 
the true neutrino luminosity. A detailed numerical calculation (without considering the 
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uncertain gravitational redshift, however) shows that the true average neutrino luminosity 
can be much higher if there is rapid rotation and/or magnetic fields holding up the collapse 
(Shi & Fuller 1998). The neutrino energy loss rate scales steeply as Tg, and the temperature 
distribution in the homologously collapsing core (an index n = 3 polytrope) follows the 
Lane-Emden function and so peaks at the center. Compensating this feature will be the Rf, 
dependence of the above energy deposition rate Qup- Therefore, we will approximate 
the entire neutrino emissivity of the core as arising from the edge of the core {Ry ~ Ts), 
and then take L^p as the characteristic neutrino luminosity from eq. (^. (Note that this 
equation is appropriate in the case where M|^*-^ ^0.1 and neutrinos diffuse from the core. In 
this case, the central temperature is irrelevant, though we may get luminosities comparable 
to the free streaming case because the core will have lower mass and, hence, a generally 
higher temperature scale.) 

The expected near-thermal spectrum of the neutrino emission implies 
{El)/{E^) = {El)/{Ep) ^ Q{M^^)-^/^MeV (Shi & Fuller 1998). Therefore, the 
neutrino energy deposition rate per unit volume will be roughly 

Q,p(r) ~ 4 X 10'' {Ml'^)-''-^{rJr)^eTgcm-h-\ (4) 

The total energy deposited into the fireball above a radius r is 

^f.b.(r) = ts / 47rr2Q,p(r)dr ~ 2.5 x 10^^ {M^'^)-''-^{r,/r)'' erg, (5) 

J r 

which is tremendous. The fireball will undoubtedly lose some of this energy to thermal 
neutrino emission. But, once the pair density is high enough for this, neutrino/electron 
scattering should deposit even more energy. If Af|^*-^ = 0.5, the energy deposited in the 
fireball will be ~ 10^^ erg at a radius r ~ Srg ~ 10^^ cm. This is the total observed energy 
in a GRB assuming a Air solid angle and a redshift 2 ~ 3. 

A successful model of GRBs must avoid excessive baryon loading so that a Lorentz 
factor of r 10^ can be achieved for the baryons accelerated by the fireball. This suggests 
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that the region at several Schwazschild radii from the supermassive star core should have 
extremely low baryon density. This may be satisfied if the whole star collapses homologously 
into a black hole, and/or substantial rotation causes the star to collapse in a flattened 
geometry with very little material in the polar directions (an extreme case of this geometry 
was discussed in Bardeen & Wagoner 1969). The homologous collapse of the entire star 
could be engineered only if the star has substantial centrifugal support from rotation 
and/or if there is significant magnetic pressure (but not so much that an explosion results). 
Therefore, rotation could be a crucial factor in this picture. Rotation will also result in a 
longer collapse timescale, and mildly beamed 7-ray emission. A high angular momentum 
collapse may therefore be challenged in generating GRBs with durations ^ 1 second. 

Another means to avoid excessive baryon loading may be possible in the collapse 
of a dense star cluster. In this case the whole star cluster can collapse on the General 
Relativistic instability (Shapiro & Teukolsky 1985) and collisions of ~ M© stars could 
provide the neutrino "engine" that powers fireballs. During the collapse, the central stars 
will have relativistic speeds and the typical entropy per baryon produced in zero impact 
parameter collisions of these will be 5 ~ lO^rV2(^^/5.5)i/4(M0/M,)i/^(K/V^0)^/^ with 
Tg ~ 1, conditions commensurate with those required for hydrostatic supermassive stars 
(S Ri lO^(M'°'VlO^M0)i/2) (jjgj-g r ~ 1 is an appropriate Lorentz factor, and VJVq is the 
ratio of the stellar collision interaction volume to the solar volume.) In fact, most collisions 
will not be "head-ons," but rather will involve the tenuous outer layers of the stars. The 
lower densities involved will translate into larger entropies (effectively, (V^/VqY^^ could be 
considerably larger), possibly large enough (S ~ 10'') to provide a pair fireball directly. In 
the collapse, space between moving stars may provide baryon-free "lanes" , and the stellar 
collisions themselves may cause the neutrino emission to be "spiky" (the overall emission 
profile, however, should nevertheless follow the free fall collapse profile indicated above for 
supermassive stars). Both processes are stochastic, possibly contributing to the "spiky" 
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time structure of the GRBs. This direct collapse of relativistic star clusters and the collapse 
of supermassive stars may well represent two extremes on a continuum of supermassive 
object collapse. 



3. Event Rate and Peak Flux Distribution 

The rate of supermassive object collapses should be able to match the observed rate of 
GRB events (several per day) if a substantial fraction of the burst events are to come from 
this source. Assuming that supermassive objects all form and collapse at a redshift the 
rate of these collapses as observed at the present epoch is 

3drpbF(l+^ 
Mio M-t ' 

where r is the Friedman-Robertson- Walker comoving coordinate distance of the objects, 
is the scale factor of the universe at the epoch corresponding to a redshift z (with ag = 1), 
to is the age of the universe, pb ~ 2 x 10~^^ Ob^^ gcm"^ fa 5 x 10~^^gcm~^ (Tytler & 
Buries 1997) is the baryon density of the universe today, h is the Hubble parameter in 100 
kms~^ Mpc~^, and F is the fraction of baryons that were incorporated in supermassive 
objects. For 2; ~ 3 we will have r ~ 3000/i~^ Mpc. The collapse rate is therefore 

0.15F [hlf^Y^ sec-^ ~ lO^F [Mf^y^ day-\ (7) 

With F ~ 0.1%, i.e., with 0.1% of all baryons having been incorporated into supermassive 
objects of ~ 10, we should observe (assuming a 100% detection efficiency) one collapse 
per day if they emitted 7-rays into a 47r sohd angle. This would constitute a substantial 
fraction of the observed rate of GRB events. The baryon fraction F = 0.1% in ~ IO^Mq 
black holes implies a (cumulative) density of 7/;,^ such supermassive black holes formed in 1 
Mpc^. This GRB rate is about two orders of magnitude lower than 24 Gpc~^yr~^, the rate 
required if GRBs originate from source populations that do not evolve over time (Fenimore 



- 11 - 



and Bloom 1995). This shortfall in rate results because we have assumed that all GRBs 
are high redshift collapse events and are therefore seen from a larger volume. In addition, 
the rate of supermassive object collapses required in our GRB model does not depend on 
the mass scale of the collapsing objects, although the fraction F scales linearly with M™*. 
Observations show that almost all galaxies that have been examined appropriately seem to 
have supermassive black holes in their centers (van den Marel et al. 1997). It is therefore 
intriguing to estimate the rate of supermassive object collapses required by our GRB model 
on a per galaxy basis. If such supermassive object collapses occurred only in normal ~ L* 
galaxies, the rate needed is about 350/i~^ per galaxy. However, this number of events 
per galaxy is much lower, perhaps ^ 10/i~^ per galaxy (based on, for example, the galaxy 
number densities of Zucca et al. 1997), if dwarf galaxies harbor supermassive objects as 
well. Therefore, it may be conceivable that these supermassive object collapse events are 
tied to the supermassive black holes at the centers of galaxies, if such supermassive black 
holes occur in every galaxy-scale object. Such an association of supermassive objects and 
galaxy-scale objects may also be born out by considering Lyman limit systems and damped 
Lyman-a systems, which are associated with galactic halos and disks at high redshifts. 
Using a column density A^hi distribution per unit column density per unit absorption 
distance of lO^^'^A^^f^^i^-^^ (Storrie-Lombardi, Irwin & McMahon 1996), we find that the 
rate of supermassive object collapse matches that of GRBs if every Lyman-a system with 
^Hi ^ 10^^ cm~^ harbors a supermassive object. 

If all GRBs are from 2; ;^ 1 then the 7-ray burst peak flux distribution (log A'"-log P) 
will be very different from models with a homogeneously distributed population of GRBs. 
The observed log A^-logP distribution is a power law with index = — 1.5 which has a break 
at the faint end (Fenimore et al. 1993). This would be consistent with homogeneously 
distributed cosmological sources with a cut-off at high redshifts, unless the peak flux of 
GRBs, P, cannot be regarded as a standard candle. But since the logA'"-logP distribution 
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is a convolution of the peak flux and spatial distribution, there is no guarantee that the 
observed power law requires a homogeneous distribution of sources. For our model, in which 
supermassive object collapses most likely occur at cosmological distances with z ^ 1, we 
can always invoke variances in the peak flux of GRBs, and/or an evolution of supermassive 
object co-moving number densities, or invoke another population of GRBs, to fit the 
observed 7-ray burst peak flux distribution. It is worth noting that even in existing stellar 
remnant-based models, the sources tend to be more abundant at ^ ;^ 1, because the star 
formation rate was higher then (Totani 1997). 

4. Conclusion 

The formation of the supermassive black holes inferred in AGNs, quasars and many 
galaxies may well involve the collapse of relativistic star clusters which form intermediate 
phase supermassive stars. We point out here that collapses of these supermassive objects 
will be accompanied by prodigious thermal neutrino emission which could transport a 
fraction of the gravitational binding energy of these objects to a region(s) where the 
baryon loading is low, thus creating "clean" fireballs that generate 7-ray bursts. The 
major advantage of this model is a huge energy release, and just such an energy scale is 
required by recent observations of high redshift bursts. We have shown that the collapse 
timescale and expected collapse event rates are consistent with 7-ray burst parameters. The 
principal weakness of our model is the baryon loading problem. We have outlined possible 
ways to circumvent this problem by appealing to high angular momentum and flattened 
collapses, and by appealing to the stochastic nature of stellar collision-induced supermassive 
star/black hole build-up in the collapse of relativistic star clusters. 

We thank David Band and Edward Fenimore for valuable suggestions. This work is 
supported by NASA grant NAG5-3062 and NSF grant PHY95-03384 at UCSD. 
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Figure Caption 



Figure 1. A flow chart for tlie collapse of supermassive objects. 
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